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A planar Mott insulator with easy plane Neel order can be mapped onto a Gutzwiller-projected 
spin Hall state. The in-plane substitution of 7aV?^ and Li^ for Cu^^ brings zero modes around the 
spin vacancy. The Gutzwiller projection precludes single particle excitations, however, preserves the 
zero modes which have a local spin moment for Zn^^ and a local charged hole for Li"*", respectively. 
While the local spin moment for Zn^^ is screened by the long-range spin correlations, the active 
charge degree of freedom for Li""" impurity twists the spin background with long ranged disturbances. 
This proposal explains the very different roles of the Zn^+ and Li"'" impurities on the magnetic 
perturbations in the Cu02 plane in the undoped cuprate compounds. 



The parent compound of high Tc superconductors, e.g. 
La2Cu04, is antiferromagnetic (AF) insulator in which 
the long range AF order has the easy plane anisotropy 
in the Cu02 plane due to the spin-orbit couplings |T]- 
[3j. It is generally believed that the high T^. supercon- 
ductivity arises out of the coherent charge behavior of 
the dopant holes in the short range AF ordered spin 
backgroundjH [5]. For La2Cu04, a very small concen- 
tration of holes 2%), introduced by Sr or Ba, suffices 
to suppress long range AF order and then the super- 
conductivity appears immediately. The insulator- metal 
transition is accompanied by different competing mag- 
netic and charge orders making it difficult to straighten 
out the magnetic disturbance from the dopant. The in- 
plane substitution of Zn^"*" and Li+ for Cu^"*" provides 
the opportunity to study magnetic disturbances from the 
dopant without competing charge orders [5]. The Zn^''" 
substitution brings the spin vacancy to the plane and 
disturbs the spin magnetism only with short ranged per- 
turbations. The long range Neel order for the Zn^+ sub- 
stitution survives until the concentration close to the 
site dilution percolating threshold [7|. Comparing with 
Zn^+ substitution, Li"*" brings an extra hole with it. The 
Li"*" substitution shows a rapid suppression of long-range 
AF order like Sr but without a transition to a conduct- 
ing state[51 [3]. The commensurate magnetic correla- 
tions in in-plane Li"*" doped La2Cu04 are very different 
from the incommensurate ones in Sr^+ or Ba^+ doped 
compounds |10j. The hole charges introduced by the Li"*" 
are well localized without the competing charge stripe 
order, however, disturb the magnetism with long range 
perturbations. 

In this letter, we propose a constructive approach to 
investigate the different roles of the non-magnetic Zn^"*" 
and Li"*" impurities on the magnetic perturbations in the 
Cu02 plane in the undoped cuprate compounds. We im- 
plement the Gutzwiller-projected spin Hall (GSH) state 
to model the planar Mott insulator with long range easy 
plane Neel order. The specified spin Hall (SH) state is 
invariant under the time reversal and ^-axis spin Sz ro- 



tation symmetries. The Gutzwiller projection freezes the 
charge degree of freedom and leads to long range spin 
correlations with the XY AF spin order jlTJ [12]. The 
non-magnetic impurity substitution brings zero modes 
for the unprojected SH state which have a local spin mo- 
ment or local charge [T5]. Due to the spin Sz rotation 
symmetry, the spin Chern number is well-defined IT5] . 
It is found that the topology described by the spin Chern 
number is invariant under the Gutzwiller projection. In 
this letter, GSH is taken as the candidate state for the 
spin system with easy AF order in the Cu02 plane. The 
in-plane Zn^"*" and Li+ substitution states are taken as 
the Gutzwiller-projected non-magnetic impurity doped 
SH states with a local spin moment for Zn^+ and a local 
charged hole for Li"*", respectively. The Gutzwiller projec- 
tion pushes the single particle excitation into the infinite 
energy region, however, preserves the local spin moment 
and charged hole represented by zero modes. The spin 
Chern number, = 2 for all concerned states, brings 
the mutual Chern-Simons term that captures the rela- 
tions between the charge and magnetic degrees of free- 
dom. While the local spin moment for Zn^"*" is screened 
by the long-range spin correlations, the active charge de- 
gree of freedom for Li"*" impurity generates the vortex tex- 
tures which destroy the long range magnetic order with 
the small critical doping around ^ 0.025. 

We start from the unprojected SH state of the tight 
binding model at half filling 

-ffsH = ^ i y fia fj<y, ( 1 ) 

ij 

with the non-zero nearest neighbor (n.n.) and next n.n. 
(n.n.n.) hopping parameters on the square lattice (Fig. 
[l} . We double the unit cell by inserting uniform ±7r flux. 
The hopping parameters, e.g. in the plaquette □1234 in 
Fig. [1] are chosen as the complex numbers with the fluxes 
as 

with — — ^^S.itijtjk^'ki) ^^'^ cr — ±1 in formula ex- 

pressions while a =t / i in the index and subscript. The 
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FIG. 1. (Color online) The square lattice without and with in- 
plane impurity substitution for (a) and (b), respectively. The 
blue square marker in (b) is the vacancy prohibiting hopping 
process on the blue dash-dot links. 



in-plane substitution of non-magnetic impurity prohibits 
the hopping process onto the spin vacancy site which is 
marked as the blue square in Fig. [l] (b), i.e. all the 
hopping amplitudes on the blue dash-dot bonds are zero. 

We impose the general twisted boundary phase on the 
finite size system, f^+L^a = ftae'^'ji+Lya = /i^e'''«. 
Under the twisted boundary conditions, the tight bind- 
ing model is easily diagonalized. In this letter, the nu- 
merical calculations are carried on the square lattice on 
the torus with the size — Ly — 8 and the hopping am- 
plitudes \t\ = 1 and |<i| = 0.2|i| for the n.n. and n.n.n. 
bonds, respectively. The energy revolutions with twisted 
boundary phase {ky = fixed) are shown in Fig. [2] 
(a) and (b) for the undoped and the in-plane impurity 
substitution SH model, respectively. For any boundary 
twist kx,ky, we can construct the undoped SH state as 

|*Sh) = |*t) a I*;) with 1^-,) = \el) A • • • A I J/'). Here 
wedge product A denotes the fermionic antisymmetry for 
the exchanges of two electrons. {|e^)} (* = 1: ' ' ' 7 ^/2) 
are the negative energy levels {N = x Ly). The va- 
cancy brings the zero modes that have four different va- 
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cancy doped SH states [T^ 

|*^^H';e) = |*nA|*r)A|Ot>A|Oj, (3) 

with l*™'^) ^ \el) A ••• A |e^/^"^). The zero modes 
{|0^), |0^)} bring the charge quantum number Q — ±1 
{Sz — 0) for the local charges in l^'g^ ; o/e) and the spin 
quantum number = ±1/2 (Q = 0) for the local spin 
moments in l^'gnift / -IJ-)j respectively [13]. 

Since Sz is conserved, it is convenient to define the spin 
Chern number Cg O [15] 

Cs = -^ I dkx [ dkyQ{kx,ky), (4) 

with the spin Berry curvarture i~l{kx,ky) — Vk x 
(^'jiVkl^') under the twisted boundary phase k^. — 
—k^ = kx,k^ = = ky. All the states have the to- 
tal spin Chern number = 2. So the topology of the 
SH state is robust against the vacancy. In Fig. [s] (a), we 
shown the fc-dependent spin Berry curvatures rt{kx,ky) 
for the state, e.g. |\I'sh;o)- It is useful to characterize 
the topology of the SH state in terms of the response to 
two different gauge fields a"^ and 

^ = / drJ2fU9rS^, + tr,e'('^"~^+'^')-"--- )/.<., (5) 

which leads to the effect Lagrangian for the gauge fields 

1 C 
g Ztt 

The first is the Maxwell term and the second is the mu- 
tual Chern-Simons term with the coefficient ^ — -. As 
emphasized in Refs. |ll|and[T^ a fiux quantum of a* car- 
ries the spin quantum number Sz — ^ and the magnetic 
field of a" is the Sz density. Meanwhile, a flux quantum 
of a'^ will create (or annihilate) one pair of spin-up and 
spin-down electrons. It carries the charge quantum num- 
ber Q — 2 and the magnetic field of a'^ is the density of 
the charged pairs. 




FIG. 2. Energy revolutions with twist boundary phase ky 
{kx — fixed) for the undoped (a) and impurity doped (b) 
spin Hall states, respectively. 



FIG. 3. (Color online) Spin berry curvature ^l{kx,ky) vs the 
twist kx and ky for (a) unprojected j^§ff;o), (b) projected 
i*GSH;o), respectively. 
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We obtain the GSH state by projecting out the double 
occupied state 

I^gsh) =7'g|%h), rG^Y[{l-n,^n,^). (7) 

i 

Here Vg is the fuUy Gutzwiller projection operator. The 
spin Chern number is unchanged under the Gutzwiller 
projection. The role of the projection is to freeze the 
charge degree of freedom of the electrons while the 
physics of the spin remains. Under Gutzwiller projec- 
tion, is removed in |^gsh) and the spin fluctuation 
is controlled by the Maxwell term a'^. Due to the 2:- axis 
spin Sz symmetry, the instanton effect is not important 
and the gauge field a" is always in the Coulomb phase 
and gives rise to a linear dispersed gapless photon mode 
corresponding to the gapless spin wave excitation. The 
gapless spin density fluctuations imply that e'^^^' spin ro- 
tation is spontaneously broken. So I^I/gsh) has the long- 
rang XY spin order which has the AF pattern directly 
confirmed by the numerical calculation [iT]. 

In this letter, we will take Iv^gsh) as the approxi- 
mated candidate state for the spins in the Cu02 for 
the cuprates. The in-plane substitution state is writ- 
ten as Wjg^j.pG I *Sh) with the vacancy inserting operator 
VFjaj.. In this paper, we assume the permutation, i.e. 
[Vg, W^Jac] = 0; so that we have the in-plane substitution 
state 

|vi/-h;*)^7'gI^?h^;*), (8) 

with * G {o, fl", -IJ-}. The Gutzwiller projection pushes 
the single particle excitation into infinite energy re- 
gion, however, preserves the local spin moment and hole 
in l^gy ; -ff, o). So the permutation assumption, i.e. 
[T'cjW^Jac] — O7 is quite reasonable. We will study the 
in-plane substitution states I^G^gHitr) and I^'gshS-IJ-) to 
investigation the magnetic disturbances for Zn^+ and 
l^(3SH' ^o'' Li^i respectively. The projected states have 
the same spin Chern number as unprojected ones, i.e. 
Cs = 2. In Fig. [3] (b), we show the fc-dependent spin 
Berry curvature for the projected state, e.g. I^'gsh! o). 

The Zn^+ in-plane substitution states, i.e. l^'cgH^ft) 
and 1 4* ; -II) I have no charge degree of freedom af- 
ter the Gutzwiller projection, i.e. Q = 0. The zero 




FIG. 4. The spin Sz distribution for the in-plane substitution 
state |*2fsH;^)- 



modes brings a local moment with the spin Sz distri- 
bution shown in Fig. [4] for l^'^sH^'fr)- 1^ is localized 
around the spin vacancy in the staggered pattern arising 
from the Gutzwiller projection. In the thermodynamical 
limit, the spin fluctuation is described by the dynamics 
of the gauge fields a'* in the Maxwell term 

■^MW — -^{^tj.iy\d^al)^, (9) 

with the gapless excitation implying the long range XY 
spin order. In the continuum limit, we have the static 
equations of motion 

V^a^ = 0, V X a^v) = z2^(5f,,(r) + ^^^.(r)), (10) 

with Qq = 0. 5'io(.(r) is the spin Sz distribution as shown 
in Fig. [4] and can be approximated as ^'[^^.(r) = ^S{r) in 
the continuum limit. Si^(r) is the induced spin configura- 
tion that screens the local moment. When S^ excitations 
are gapless, we can obtain the full screening for the lo- 
cal moment resembling the Kondo screening for the local 
moment in the metal. The local moment will survive 
provided that the spins are short range correlated and 
the screening is not complete. This may be the reason 
for the local moment associated with the dopant Zn^+ 
in the pseudogap phase with short range spin order in 
the underdoped cuprate compounds. Regardless of the 
screening, the dopant Zn^"*" brings only the short range 
magnetic disturbance and doesn't change the XY AF 
spin order in I^gsh;!!" / -U-)- It behaves as the site dilu- 
tion percolation. 

Comparing with Zn^+ , the hole in the in-plane substi- 
tution state of Li"*", i.e. I^'cgn! o), replaces the local mo- 
ment in |^'™gH;'fr / -IJ-)- The distributions of the hole and 
spin for I^'qsh; o) are shown in Fig. [5](a) and (b), respec- 
tively. Due to the correlated effect of the Gutzwiller pro- 
jection, the distributions vary from S^{r) for |^'™sH;'fl") 
shown in Fig. [ijto nh{r) for |^'gsh;o) Fig- [5](a). It 
is intriguing that the charge degree of freedom for Li+ 
activates the mutual Chern-Simons term 

-S^MCS — —^uv\o!^iidya\. (11) 

In the continuum limit, we have the equations of motion 

V^ag =7i,,(r), a^=0, (12) 




FIG. 5. (a) The hole Uh and (b) spin 5*^ distribution for 

I*gsh;o). 
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where n;i(r) = :^eo^i/9^a^ is the charge distribution of 
the hole in | ^'™gH i '-') 

Fig. [5] (a). There is no induced 
charge that screens n/i(r). Due to the mutual Chern- 
Simons term, the charged hole generates the tt and — tt 
fluxes in fi-^ and fi^, respectively. Therefore, Li+ gen- 
erates a 27r vortex in the XY spin order 5+ = fi^fii- 
Considering the compactness, the gauge field should 
be written as 



(13) 



with the integer field A/"^ {n e Z) that keeps track of 
the compactness. When n ~ —1, the integer field A/'^ 
can convert the 27r vortex in XY spin order into —2tt 
anti-vortex. 



From the mutual Chern-Simons term (ff ), we know 



that Zn^+ does little damage on the XY AF spin or- 
der while Li"'' generates the destructive vortex structure. 
Below we shall confirm the results by numerical compu- 
tation directly. For the finite system, we have the expres- 
sion for the AF spin order pH IT^ 



Map = (-l)'"'(^'GSH|^+|*GSH.(-l)-flux) 



(14) 



with the projected minus 1-flux state |^GSH,(-i)-flux)- 
(— 1)''' is the AF prefactor. Without in-plane substi- 
tution, the spin configuration is in the AF Neel patter, 
Maf ~ const. pTj. For Zn^+ substitution, the AF spin 
order for l^'csH^tl") is given as 

A^If = (-1)""-(*gsh;^ |5+l^lsH,(-i)-flu.), (15) 
with the projected minus f-flux state 



l*^SH,-i)-flux) = -Pal^iTn A IvPlX"^), (f6) 



-i)-flij 



with l^f^zlf = lei) A • • • A lef/'"') and |* 



(-l)-flux> 



|e|) A • • • A |e^^^)- |e^) is the negative energy level in the 
presence of the uniform flux integrated to a minus flux 




,^ 

. /' 




FIG. 6. XY AF spin configuration MAF(ri) for (a). Zn^+ 
substitution state |^GSHi1^)i (b). Li^ substitution state 
I^J'S'sHi o)j respectively. The AF prefactor (— 1)"^' is already 
included in the definition of MAF(ri), e.g. in Eq. (141. In 



(b), the imposed periodic boundary condition on the torus 
generates the anti-vortex paired with the vortex generated by 
the Li+ impurity 



quantum. Similarly, the XY AF order for Li"*" substitu- 
tion state I^'osh! o) is given as 

= (-l)'"-(*GSH;o|S.+ l*GSH,(-l)-flux), (17) 

with the projected minus 1-flux impurity doped SH state 
l^-^H.t-D-flux) = ^Gl^^il"'"'') A (18) 



.(-1)- 

with l^l-.i^:""'^) = lei) A 



A e 



|e|)A 



A|e/ \ 



N/2-2 
t 



) and Y^'f^'-) = 



In Fig. [H we show the numerical results of XY AF spin 
order for Zn^"*" and Li"*" substitution states in (a) and (b) , 
respectively. As expected, I^qshSI)") h^-s the untwisted 
XY AF spin order while | 5'™sh ; o) has the twisted spin 
order with the vortex structure. It should be noted that 
the imposed periodic boundary condition on the torus 
generates the anti-vortex paired with the vortex gener- 
ated by the Li+ impurity. 

The vortex structure generated by Li+ substitution can 
destroy the long range spin order due to the Kosterlitz- 
Thouless transition [TS]. In the presence of the vortex- 
antivortex pairs, the renormalized spin stiffness at tem- 
perature T is given as [T7Hl9] 



T 



|-4-V(f) 



dr 



3-2irps/T 



(19) 



with ps is the bare stiffness without vortex and = 
Vhi + Vth the total fugacity of vortex pairs where 
for the pairs generated by Li"'" and y'^^ for the thermally 
generated pairs, tq is the size of the vortex core. At 
the critical temperature Tc, we have = 2Tc/tt. The 
rigidity of the spin order can only sustain the amount of 
the vortex-antivortex pairs with the fugacityfUl [19] 



yl + yk-^{i~—)'- 



(20) 



At the critical concentration of Li"'", Tc = and y^j^ — 
such that y£; = and we obtain the critical concentra- 
tion 



1 

4^ 



= 0.025, 



(21) 



with the estimation tq = 2a. 

During our construction, all the unprojected states are 
protected by the finite gaps. The spin Chern numbers 
of these states are well-defined and remain unchanged 
after projection. The key point in our construction is 
the nonzero spin Chern number which leads to the mu- 
tual Chern-Simons term when the charge degree of free- 
dom is activated. The mutual Chern-Simons structure 
in the doped antiferromagnet is extensively studied for 
the cuprate by Weng's group starting from the bosonic 
resonant valence state[TlH22]- The spin and charge tex- 
ture for the substitution of Li+ in the Cu02 plane was 
also studied in Ref. 1231 where the skyrmion structure was 
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proposed. The spin texture induced by the delocaUzed 
hole was also studied in the the eary literatures of high 
Tc cuprateslMlES]. 



In conclusion, we present a constructive approach to 
investigate the magnetic disturbance for the in-plane sub- 
stitutions for Zn^+ and Li''" for Cu^+ in the Cu02 planes 
in the undoped cuprate compound. In the substitution 
states, Zn^+ and Li"*" impurities can be represented as va- 
cancies introducing a zero mode, which has a local spin 
moment for Zn-^+ and a charged hole for Li+, respec- |~ 
tively. While the local spin moment for Zn^"*" is screened 
by the long-range AF spin correlations, the active charge 
degree of freedom for Li+ impurity generates the vortex 
in the spin AF background and suffices to suppress the 
spin order at the critical concentration Xc ^ 0.025. 



During the work, from the issue raising to the 
manuscript preparing, many helpful suggestions and 
comments have been given by T.M. Rice. JW Mei also 
thanks M. Sigrist and Z.Y. Wcng for useful discussions. 
The suggestion of numerical methods from Lei Wang is 
also useful. The work is supported by Swiss National- 
Fonds. 



[1] 



C. J. Peters, R. J. Birgeneau, M. A. Kastner, 
H. Yoshizawa, Y. Endoh, J. Tranquada, G. Shirane, 
Y. Hid aka, M. Oda, M. Suzuki, and T. Murakami, jPhys] 
Rev. B 37, 9761 (1988) 

D. Coffey, T. M. Rice, and F. C. Zhang, |Phys. Rev. B 
44, 10112 (1991) 

T. Yildirim, A. B. Harris, A. Aharony, and O. Entin- 
Wohlman, Phys. Rev. B 52, 10239 (1995), 
P. W. Anderson, Science 235, 1196 (1987 



_P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys, 
78, 17 (2006); 



[6: 

[7 
[9 

[lo: 



[11 

[12 
[13 

[14 

[is: 
[16: 



H. Alloul, J. Bobroff, M. Gabay, and P. J. Hirschfeld, 
Rev. Mod. Phys. 81, 45 (2009) 

O. P. Vajk, P. K. Mang, M. Greven, P. M. Gehring, and 



J. W. Lynn, Science 295, 1691 (2002)[ 
A. I. Rykov, H. Yasuoka, and Y. Ueda, |Physica C: Su-] 
Ijerconductivity 247, 327 (1995) 

J. L. Sarrao, D. P. Young, Z. Fisk, E. G. Moshopoulou, 
J. D. Thompson, B. C. Chakoumakos, and S. E. Nagier, 
|Phys. Rev. B 54 , 12014 (1996) 

W. Bao, R. J. McQueeney, R. Heffner, J. L. Sarrao, 
P. Dai, and J. L. Zarestky, |Phys. Rev. Lett. 84, 3978| 
(2000) 

Y. Ran, A. Vishwanath, and D.-H. Lee, ArXiv e-prints 
(2008), arXiv:0806. 2321 [cond-mat.str-el] 
Y. Ran, W.-H. Ko, P. A. Lee , andX.-G. Wen, [ Plyl] 
Rev. Le tt. 102, 047205 (2009)| 

J. He, Y.-X. Zhu, Y.-J. Wu, L.-F. Liu, Y. Liang, and 
S.-P. Kou, ArXiv e-prints (2012), arXiv: 1210.0266 

D. N. Sheng, Z. Y. Weng, L. Sheng, and F. D. M. Hal- 
dane, Phys. Rev. Lett. 97, 036808 (2006)l 

E. Prodan, Phys. Rev. B 80, 125327 (2009)] 

J. M. Kosterlitz and D. J. Thouless, Journal of Physi"cs| 
C: Solid State Physics 6, 1181 (1973) 



P. Chaikin and T. Lubensky, Principles of Condensed 
Matter Physics (Cambridge University Press, 2000). 
C. Timm and K. H. Bennemann, ,Phys. Rev. Lett. 84, 
4994 (2000) 

_J^W. Mei and Z. Y. Weng, Phys. Rev. B 81, 014507 

(2010) 

S.-P. Kou, X.-L. Qi, and Z.-Y. Weng, |Phys. Rev. B 71, 
235102 (2005) 



P. Ye, C.-S. Tian, X.-L. Qi, and Z.-Y. Weng, |Phys. Ri^ 
Lett. 106, 147002 (2011), 

P. Ye, C.-S. T ian, X.-L. Qi, and Z.-Y. Weng, [Nuclear 
Physics B 854'^, 815 (2012) 

S. Haas, F.-C. Zhang, F. Mila, and T. M. Rice, |Phys. 
Rev. Lett. 77, 3021 (1996) 



B. I. Shraiman and E. D. Siggia, Phys. Rev. Lett. 61, 
467 (1988) 



B. I. Shraiman and E. D. Siggia, Phys. Rev. B 42, 2485 



(1990)1 



